Detecting entanglement in two mode squeezed states by particle counting 
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We present an entanglement criterion for two mode squeezed states which relies on particle count- 
ing only. The proposed inequality is optimal for the state under consideration and robust against 
particle losses up to 2/3. As it does not involve measurements of quadratures - which is typically 
very challenging for atomic modes - it renders the detection of atomic many-particle entanglement 
feasible in many different settings. Moreover it bridges the gap between entanglement verification 
for a qubit and criteria for continuous variables measured by homodyne detection. We illustrate its 
application in the context of superradiant light scattering from Bose Einstein condensates by con- 
sidering the creation of entanglement between atoms and light as well as between two condensates 
in different momentum states. The latter scheme takes advantage of leaving the Gaussian realm 
and features probabilistic entanglement distillation. 

PACS numbers: 



I. INTRODUCTION 

Entanglement is a true quantum feature. The study 
of this peculiarity of physics does not only hold the 
promise to acquire a deeper understanding of Nature, 
but also paves the way towards auspicious applications 
of quantum information science such as unconditionally 
secure communication, ultraprecise measurements, quan- 
tum computing and quantum simulation. Therefore, the 
quest for entanglement or inseparability criteria is a vig- 
orous field of research [l[ . 

In the most basic case, entanglement is shared between 
two parties holding a single particle each and bipartite 
entanglement of single pairs is well understood. Then 
again, entanglement between parties holding a huge num- 
ber of particles was studied for Gaussian states with great 
success 0]. Here, the natural question arises, how entan- 
glement can be verified in the intermediate regime and 
how the two well-studied cases cases of single pairs and 
Gaussian modes can be linked by an inseparability cri- 
terion for an arbitrary number of particles, which is not 
restricted to Gaussian states. Starting from this motiva- 
tion, we bridge the gap between inseparability criteria for 
a qubit and continuous variables for entanglement in two 
mode squeezed states. Moreover, from a practical point 
of view, the new entanglement condition provides experi- 
mental feasibility and applicability in numerous settings. 

More specifically, entanglement in two mode squeezed 
states can, in principle, be detected by means of a Gaus- 
sian inseparability criterion 0], which requires the mea- 
surement of variances in canonical quadratures. This 
can be conveniently performed for light modes, as first 
demonstrated in (4J, as well as for multi-atom collec- 
tive spin modes, as shown in However, in many 
cases involving multi-atom entanglement homodyne mea- 
surement of atomic canonical variables which require an 
atomic "local oscillator" is not feasible. 

This problem can be overcome by means of the prac- 
tical entanglement condition put forward in this work. 



The proposed inequality requires only particle number 
measurements, rather than measurement of quadratures 
and can be used to detect TV-particle entangled states of 
the form \TMSS) <g) \TMSS), which is attractive in view 
of many recent experiments, that offer the potential to 
generate this type of entanglement. 

For example, superradiant scattering Q of laser light 
from a Bose-Einstein condensate was observed recently 
Q. Superradiant scattering leads to highly directional 
emission of light from the atomic sample. This striking 
effect attracted considerable interest and it was shown 
theoretically that two mode squeezed states can be gen- 
erated in this context This is particularly inter- 
esting, since it represents the interspecies atom-light ana- 
logue of photonic twin-beams generated in optical para- 
metric down conversion, which plays an essential role in 
many applications of quantum optics and quantum in- 
formation theory. Despite the fact that the dynamics of 
the process and the resulting non-separable state are well 
understood and the system is known to be a very promis- 
ing candidate for entanglement generation d, [l(| [ll[ , 
entanglement could not be verified in this system due to 
the absence of a suitable inseparability criterion. Thus, 
nonclassical correlations have been studied, but since the 
quantum states produced in an experiment cannot be as- 
sumed to be pure, correlations do not imply the presence 
of entanglement. Other examples can be found in many 
different setups, as for example in entanglement produc- 
tion in spin exchange collisions in Bose-Einstein conden- 
sates [l2| , the free electron laser [l3[ , the creation of two 
mode squeezed states by dissociation [l4| and entangle- 
ment production in colliding Bose-Einstein condensates 
[l5| or four-wave mixing in matter waves in an optical 
lattice 16]. While correlations could be observed, multi- 
particle entanglement has not yet been detected in this 
context. This gap can be closed by means of the entan- 
glement criterion put forward in this work. It belongs to 
the class of entanglement conditions derivable from the 
partial transposition criterion [ljj and as the Gaussian 
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criteri on [3| , it can be seen as as a local uncertainty re- 
lation with the difference that it involves a bound, 
which is given by expectation values of operators, rather 
than uncertainty limits represented by fixed numbers. 
The inequality is optimal for the state under consider- 
ation. Moreover it is robust against the sources of noise 
to be expected in a realistic setup and provides a pos- 
sibility for successful detection of entanglement even for 
highly mixed states. The inseparability criterion is pre- 
sented and proven in section [TTJ In this section we also 
consider the influence of particle losses and find that for 
symmetric particle losses in Alice's and Bob's systems, 
loss of a fraction up to 2/3 of all particles can be toler- 
ated. In section IIIII we illustrate the application of the 
inseparability criterion by means of two specific exam- 
ples in the context of superradiant light scattering from 
Bose-Einstein condensates. More precisely, we consider 
the creation of entanglement between atoms and light in 
superradiant Raman scattering and describe a scheme in 
which entanglement between a moving condensate and a 
condensate at rest is created and purified. 

II. INSEPARABILITY CRITERION BASED ON 
PARTICLE NUMBER MEASUREMENTS 

As outlined above, the proposed entanglement crite- 
rion is optimal for the state \TMSS) ® \TMSS). In the 
Fock basis this quantum state is given by 

oo oc 

= (1 - A 2 ) A»N <8> A"»M, (1) 

n— m— 

where A = tanh(r) and r € C is the squeezing parameter. 
The second and fourth ket refer to Alice's system, which 
is described by two modes with creation operators 

and al- Likewise, the first and third ket refer to Bob's 
system, which is described by two modes with creation 
operators and b_. Using this notation, we change to 
the Schwinger representation and define Stokes operators 
S x , S y and S z for Alice's system and J x , J y and J z for 
Bob's. S is given by 

S x = [tia,x - n>A,y)/2, 

Sy = {UA,+45 - n A - i5 )/2, (2) 

S z = (tm,+ - tu,-)/2, 

and J is defined by analogous expressions. The 
number operators n carry subscripts for Alice's/Bob's 
system(A/ B) and for the three different bases (x/y, 
+45/ — 45, +/— ), where tia ± = a s ±a± and n B ± = &±&±- 

A. Entanglement criterion 

The main characteristic feature of the two mode 
squeezed state |*i) is the correlation of particle num- 
bers in Alice's and Bob's system. In contrast, any sep- 
arable state satisfies a lower bound for the difference in 



particle numbers in different mutually independent bases. 
More precisely, any 2x2- mode bipartite separable state, 

P = YuiViPii Pi = Pt®pf (where k > and J2iPi = 1)j 
satisfies the inequality 

(( J X ~ S X f) p + ((Jy + Sy)\ + ((J Z ~ S Z f) p 

>((n A ) p + (n B ) p )/2. (3) 

In the following, we prove this entanglement criterion. 
For any pi, the left side of inequality ((3]) equals 

(J x + Jy + J z ) Pi + (S% + Sy + S z ) Pi 
—2(J X S X — JySy + J z S z ) Pi . (4) 

(Jl + J 2 y + J 2 Z ) P , = ((n B /2) ((n B /2) + l)) Pi , where n B = 
n B ,+ +iiB-, and an analogous equality holds for the sec- 
ond term, as can be inferred from definition ([2]). Since 
Pi is assumed to be separable, the third term in expres- 
sion (j4j can be expressed as a product of two expectation 
values 

-2(J S) w = -2<J) W <S) W 

where S = ( S x -S v S z ) T . Using (J) Pi < (n B /2) Pi 

and (S) Pi < (n A /2) Pi , we obtain that expression (j4]) is 
greater than or equal to 




Since (n^/2) Pi (n B /2) Pi = (riAn B /4) Pi for product 
states, we can reexpress this equation by 

((t-t)')„ + (t + t>^<t + t> b - 

As this is true for every pi, the average ((tia + n B ) /2) p 
is a lower bound for the mixture p — YliPiPi- 

This limitation imposed on convex mixtures of product 
states can be overcome if entanglement is involved. In 
particular, 

{(Jx - ^) 2 >|*i>(*i| + ((Jy + Sy) 2 )\t, 1 )^ 1 \ 

+{(J Z - S z ) )|* 1 >(* 1 | = o, 

as the two mode squeezed state |*i) is a simultaneous 
eigenstate of (J x — S x ), (J y + S y ) and (J z — S z ) with 
common eigenvalue 0. 

B. Implications of particle losses 

In this subsection, we analyze how particle losses im- 
pair the performance of the presented entanglement cri- 
terion. The influence of particle losses is modeled by 
a beam splitter transformation. Creation operators for 
atoms and light transform according to 

al_ i-> VI - r A a\ - i^/rj v\ ± , 
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corresponding to a beamsplitter with reflectivity r A for 
Alice's system and r B for Bob's system. The quantum 
noise operators v A± , and v B± obey canonical commuta- 
tion relations for each mode and are mutually indepen- 
dent. For the target state \^x), the left side of condition 
([3]) is transformed into 

~ {(r A - r B ) 2 (An) 2 + (r A (l - r A ) + r B (l - r B )) (n»,(5) 

where (n) = sinh 2 (r) and (An) 2 = 2sinh 4 (r) is the vari- 
ance of the particle number [22j . By applying the same 
beamsplitter transformation to the right side of inequal- 
ity §5§ one obtains 

(l-r A )(n) + {l-r B )(n). 

In case of symmetric losses r A = r B = T, successful en- 
tanglement verification requires therefore r < 2/3. Limi- 
tations imposed by particle losses which are different for 
Alice's and Bob's system are more restrictive, as they 
impair directly the symmetry-property to which the cri- 
terion is tailored to. For large particle numbers, the first 
term in expression ([5]) is likely to hinder the detection 
of entanglement. This problem can be resolved by intro- 
ducing gain factors g A and g B for Alice and Bob, which 
characterize the amplification of measured signals 

(4 h-> \/g A (l ~ r A )a f ± - i^/rj v\ ± , 
b'± h-> \/g B (l - r B )b f ± - iy/r^ v B± . 
In this situation, we obtain the result 

^(g B (l-r B )~g A (l-r A )) 2 (An) 2 + 
3 

2 {9bt b {1 - r B ) + g A r A (l - r A )) (n), 
which has to be compared to 

9b(1 - r B )(n) + g A {l - r A )(n). 

g B and g A have to be optimized according to the ex- 
perimental parameters. For large particle numbers, the 
quadratic term will dominate such that g A /g B = (1 — 
— r A ) renders the problem of asymmetric losses 
and we obtain the condition (r A + r B )/2 < 2/3 for suc- 
cessful entanglement detection. Remarkably this thresh- 
old does not depend on the degree of squeezing. If the 
probabilities for particle losses are known, atomic and 
photonic signals need not to be amplified, as it is suffi- 
cient to adjust the inequality accordingly. 

III. ENTANGLEMENT IN SUPERRADIANT 
SCATTERING FROM BOSE- EINSTEIN 
CONDENSATES 

As explained above, the presented entanglement crite- 
rion can be applied in many different settings. We de- 
scribe here the verification of entanglement produced in 
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FIG. 1: Setup and atomic levels considered for the creation of 
iV-particle entanglement between atoms and light, a) Light 
propagating along z is scattered from an elongated conden- 
sate. Photons are scattered superradiantly into two endfire- 
modes, which correspond to scattering angles of and it, as 
indicated by arrows, b) Off-resonant laser light couples to 
the transition |c) — > \e), such that atoms, which are initially 
prepared in |c) are transferred to state or \b+) or |6_) via 
a Raman process and emit a photon in + or — polarization 
respectively. 



superradiant scattering of laser light from a Bose Ein- 
stein condensate. More specifically, we consider a Bose 
Einstein condensate, which is elongated along z and ex- 
cited by a laser field propagating along the same direc- 
tion. The scattering interaction is assumed to be well 
within the superradiant regime, such that light is pre- 
dominantly emitted along two endfiremodes [8|, along z 
and — z , as shown in figure [TJl. Atoms scattering pho- 
tons along ~'z acquire a momentum kick of 2hk, where 
k is the wave vector of the incoming light field. These 
atoms get spatially separated from the BEC and form a 
new moving condensate. To begin with, only the end- 
firemode consisting of photons scattered along — z and 
atoms traveling along the z direction are considered. 

A. Entanglement between atoms and light 

Atoms are assumed to possess an internal level struc- 
ture as shown in figure [TJd. As was shown in 0, the 
dynamics of the superradiant process can be described 
by a two mode squeezing Hamiltonian 

E x oc alb^+albl + E.C. , 

where the creation operators a , and a! denote the scat- 
tered light fields in plus and minus polarization, while 
b\ and bl_ are the creation operators for the respective 
atomic states. This leads to the generation of the two 
mode squeezed state Atom and photon numbers 

are correlated for each polarization and inseparability of 
the produced quantum state can be verified according to 
criterion ([3]) by identifying Alice with the light field and 
Bob with the atomic system. In the considered physical 
setting, various sources of noise may impair the verifica- 
tion of entanglement. Apart from particle losses, which 
have been discussed in the previous section, undesired 
atomic transitions can degrade the reliability of the pro- 
posed criterion, for instance when atoms are scattered 
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into states other than and |6+), while emitting pho- 
tons in + or — polarization. These processes can be 
avoided by a suitable choice of atomic levels. As an ex- 
ample we consider typical Alkali atoms used in BEC ex- 
periments, 87 Rb and 23 Na, which have nuclear spin 3/2. 
By preparing the atomic sample in the F = 1, mj? = 
ground state and inducing transitions to a manifold with 
F' = 0, mpi = 0, atoms can only be scattered back to 
the F = 1 groundstate manifold, occupying the states 
|F = l,m F = -1) = |6_) and \F = l,m F = +1) = |6+), 
while transitions to other states are forbidden due to the 
selection rule A(F) = 1. Unintentional transitions may 
also be mediated by interatomic collisions. The effect of 
transitions from \b+) or |6_) to other states is already in- 
cluded in consideration of particle losses above, but the 
creation of a pair of atoms in these two states without 
the production of the corresponding photon pair has to 
be avoided. This can be done by applying electromag- 
netic fields imposing Stark shifts on the internal levels 
such that such a transition is prohibited by energy con- 
servation. 

The measurement of Stokes operators of light required 
for the verification of the entanglement can be performed 
in a standard fashion (Tgj]. The measurement of the 
atomic collective spin projection J z can be done by count- 
ing atoms in the final states +, — with resonant absorp- 
tive imaging. The measurement of the projections J XtV 
can be performed by applying suitable radio-frequency 
7r/2 pulses to the final atomic states and then doing ab- 
sorptive imaging and atom counting. 

B. Entanglement between two condensates 

The correlations between atoms and light, that are gen- 
erated in the process described above, can be used to cre- 
ate entanglement between two condensates. To this end, 
both endfire modes are considered. The full Hamiltonian 
is given by 

H 2 (x a^jb^j + oLjblj + a i +II b\ II + a)_ n tf_ n + H.C. , 

where the subscript I refers to the backward scattered 
light and and the moving condensate, while the subscript 
II refers to forward scattered light and the condensate 
at rest. After the interaction, the backward scattered 
light field and the moving condensate are in a two mode 
squeezed state as well as the light field scattered in for- 
ward direction and the part of the condensate at rest, 
which is transferred to state \b+) or \b-), 

oo oo 

|* 2 > = (l-A) 4 5>»|n)®X>TO> 

71=0 j = 

OO oo 

<g> A m |m)|m) ®£}A'|Z)|Z), 

m=0 1=0 

where the first and second term refer to /-operators - 
atom-photons pairs in plus and minus polarization rc- 
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FIG. 2: Success probability versus entanglement which can 
be produced by means of the proposed scheme for different 
values of A. Diamonds: A = 0.7, stars: A = 0.8, squares: 
A = 0.9. The inset shows the probability of obtaining at least 
as much entanglement as was present in the input state 



spectively - while the third and fourth term refer to II- 
operators. 

The part of the resting condensate being in state |6 + ) 
or can be entangled with the moving condensate 
by means of entanglement swapping, i.e. by measuring 
EPR operators for each polarization of light modes using 
homodyne detection. However, this procedure leads to 
degradation of entanglement if non-maximally entangled 
states are involved, and a distillation step has to be per- 
formed afterwards to obtain a more useful resource state. 
It has been shown that continuous variable entanglement 
cannot be distilled using only Gaussian operations [2(| 
[23l ]. In the following, we describe therefore a scheme 
which relies on photon counting rather than Gaussian 
measurements and exhibits probabilistic entanglement 
distillation. More specifically, the moving and resting 
condensates are entangled by combining the forward and 
backward scattered light modes at a beamsplitter and 
measuring the photon numbers at both output ports. 
This has to be done for each polarization separately. We 
explain the scheme for the + polarized part of the light 
field. Analogous expressions hold for the — polarized 
part. 

By applying a balanced beamsplitter transformation 
a!™r -> (a°f + a°ffj)/y/2, a% -> (a# - a%)/y/2 , 
where aVj/a+j/ and a°$/a°?j I denote annihilation oper- 
ators of the light fields at the input and output ports of 
the beamsplitter respectively, to state (JTJ, we obtain 

OO TL 771 

|*f*> = (I — A 2 ) £ A«+«> T 2=2^££ 
— ' vn!m! ~: ~: 

71,771= U 1=0 3=0 

( < ) ( 7 ) (-W(* + i)!(rc + m-i-j)! 

1 71, i + j)\m, n + m — i — j) . 

The probability of detecting Nj photons at the first, and 
Nji photons at the second output port of the beamsplit- 
ter is Pni,Nh = (1 - A 2 )A Nl+N ". Such an event results 
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in the quantum state 



JY 



I* 



/> = ^2 k Nt,Ni t (n) \n,N- ra), 



Ni,Ni 



where N = Ni+Nn denotes the total number of detected 
photons. The coefficients kN It N I[ (n) are given by 



Wi,N 



,A n ) 



Nil 

2 F 1 (-n,-N I ,N II 



N H \ (N-n)\ 



1 



nl 

-1; 



(JV>j-n)! 



-1) 



where 2F±(a, b, c; z)/(c — 1)1 is the regularized hyper- 
geometric function. This state describes now pairs 
of atoms in in the moving condensate and at 

rest, which are referred to in the first and second 
ket respectively. For certain measurement outcomes 
this state is more entangled than | v I / i), such that 
the state can be purified by postselection. Figure 
[2] shows the success probability versus the produced 
entanglement given by the von Neumann entropy of 
the reduced density matrix of the resulting atomic 
state E{N^N H ) = £~ =0 A&, iJV „(n)ln (n)) for 

different values of A. Note that the initial state, which 
contains infinitely many terms, is truncated by the mea- 
surement process. In this way states with entanglement 



close to the maximal degree of entanglement in the 
corresponding subspace can be produced. For example, 
for Ni — 1 and Nn = 0, the maximally entangled 
state |*) = (| 1,0) + |1,0))/V2 is created. Considering 
the light field in + as well as in — polarization, the 
resulting atomic state can be detected by criterion 
([3]) after local transformation n+u i— > ./V — n + n and 
Ti-n i > N* ri-n. In this case J refers to atomic 
operators at rest and S describes the moving condensate. 

After preparation of this manuscript, Dr. C. Si- 
mon informed us that the entanglement criterion ([3]) 
presented in this work coincides with inequality (4) 
derived in 12111. 
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